Accepted for publication in ApJ 



Probing Circumnuclear Environments with the HCN(J = 3—2) 
and HCO+(J = 3-2) lines: Case of NGC 1097 

Pei-Ying Hsieh 1 ' 2 , Paul T. P. Ho 1,3 , Kotaro Kohno 4 ' 6 , Chorng-Yuan Hwang 2 , Satoki 

Matsushita 1 ' 5 
pyhsieh@asiaa. sinica.edu. tw 

1 Academia Sinica Institute of Astronomy and Astrophysics, P.O. Box 23-141, Taipei 

10617, Taiwan, R.O.C. 

Institute of Astronomy, National Central University, No. 300, Jhongda Rd., Jhongli City, 

Taoyuan County 32001, Taiwan, R.O.C. 

Harvard- Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, MA 02138, 

USA 

4 Institute of Astronomy, The University of Tokyo, 2-21-1 Osawa, Mitaka-shi, Tokyo 

181-0015 

5 Joint ALMA Office, Alonso de Cordova 3107, Vitacura 763 0355, Santiago, Chile 

6 Research Center for the Early Universe, University of Tokyo, 7-3-1 Hongo, Bunkyo, 

Tokyo 113-0033, Japan 

ABSTRACT 

We present the first interferometric HCN(J = 3-2) and HCO + (J = 3-2) maps 
in the circumnuclear region of NGC 1097, obtained with the Submillimeter Array. 
The goal is to study the characteristics of the dense gas associated with the 
starburst ring/Seyfert nucleus. With these transitions, we suppress the diffuse 
low density emission in the nuclear region. We detect and resolve the individual 
compact giant molecular cloud associations (GMAs) in the 1.4 kpc circumnuclear 
starburst ring and within the 350 pc nuclear region. The nucleus is brighter 
than the ring in both lines, and contributes to ~ 20% and ~ 30% to the total 
detected HCO+(J = 3-2) and HCN(J = 3-2) flux, within the central 1.4 kpc. 
The intensity ratios of HCN(J = 3-2)/HCO + (J = 3-2) are roughly unity in the 
GMAs of the starburst ring. However, this ratio is up to ~ 2 in the nuclear 
region. From the HCN(J = 3-2)/HCN(J = 1-0) ratio of < 0.2 in the nucleus, 
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we infer that the nuclear HCN(J = 3-2) emission might be optically thin. The 
HCO+(J = 3-2) and HCN(J = 3-2) show correlations with 12 CO(J = 3-2) and 
the 24/zm emission. The tight correlations of HCN(J = 3-2), HCO + (J = 3-2) 
and 24/zm emission in the starburst ring suggest that the dense molecular gas and 
the dust are from the same origins of star forming regions. On the other hand, 
the HCN(J = 3-2) emission of the nucleus is significantly enhanced, indicating 
mechanisms other than star formation, such as AGN activities. A self-consistent 
check of the fractional abundance enhanced by X-ray ionization chemistry of the 
nucleus is possible with our observations. 

Subject headings: Galaxies: individual (NGC 1097) - Galaxies: ISM - Galaxies: 
Seyfert - radio lines: ISM - Submillimeter: ISM - ISM : molecules 



1. INTRODUCTION 



HCN and HCO + m olecules are important dense gas tracers, w 



in external galaxies (e.g.jRickard et al. 



2006 



Henkel et alll994lP aglion e et al.lll995 



Garcia-Burillo et al 



1977 



rich are readily detectable 



Nguyen- Q-Rieu et al. 



iKohno et al. 



2003 



1989 



Jackson et al. 



Gao fc Solomon 



2004 



1993 



Riechers et al 



Montero-Castano et al 



and in our own Galactic Center (e.g., 
20091 ). With their larger dipole moments, these rotational lines have higher critical densities 
than that of CO and hence are excellent probes of the dense cores within the molecular 
clouds. However, as the fractional abundance of HCN and HCO + are usually 1000 times less 
than CO, these lines are correspondingly fainter. Nevertheless, extragalactic HCN emission 
lines conform to the well-known correlation with infrared emission over 7-8 orders of mag- 



nitude in sizescales ( Gao &: Solomon 2004 ^ Wu et al. 20051) . The HCN-IR luminosity has a 



tighter correlation than that of CO-IR (ISolomon et al.lll992l . ll997l ). This has been attributed 
to the fact that star formation occurs in dense gas rather than the diffuse gas traced by CO. 
However, in the extreme conditions of active galactic nucl ei (AGN) and sta r bursts, HCN 



can be excited by the IR radiation rather than collisi ons (lAalto et al.l 11995c IKohno et al. 



2008 


: Sakamoto et al. 


2010: Garcia-Burillo et al. 


2010), and the a 


junc 


ances of the HCN 


and HCO + can be affected by the X-ray irradiation (Tacconi et al. 


1994 


LeDD & Dalgarno 


1996 


:IUsero et al. 


2004 


Meiierink & Soaans 20051: 


Meiierink et al. 2007: 


Imanishi et al. 2007: 


Kohno et al. 2008 


: Harada et al. 2010c Garcia-Burillo et al. 


2010). Early multi-J statistical 


studies of the extragalactic HCN and HCO + were made by 


Krips et al 


• ( 


2008) to study the 



molecular gas properties in AGNs and starbursts. The abundance and excitation of these 
molecules are indeed affected by the circumnuclear environments. Here, we focus on the 
J = 3-2 transitions of HCN and HCO + , which have ~30 times larger dipole moments as 
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compared to CO, and have transition energy that are 6 times higher (~26 K) as compared 
to the J = 1-0 transitions. Our purpose is to define the utility of these lines in the study of 
the dense and warm circumnuclear regions. 

NGC 1097 is at a distance of 14.5 Mpc, and has a low luminosity-Seyfert nucleus and 
a circumnuclear starburst ring. The n ucleus is prominent in X-ray, and is associate d with 
a supermassive (~1O 8 M ) black hole (ILewis fc Eracleoud 120061 ; iNemmen et al.ll2006l ). The 
circumnuclear starburst ring shows prominent continuum emission at cm, mid-IR, and far- 
IR wavelength s, with compact kno t s, which are associated with star formation at a SFR 
of 5 M Q yr _1 (IHummel et al.lll987t ISandstrom et al.ll2010l ). The general properties of the 
circumnucle ar molecular gas of CO, HON, and HCO + were stu died with high resolution inter- 
ferometers (jKohno et al.ll200ll . 120071 ; iHsieh et al.ll2008l 1201 lh . These lower- J studies show a 
molecular ring and a nuclear concentration, which are coincident with the starburst ring and 
the Seyfert nucleus. The HCN(J = 1-0) enhancement in the nucleus rela tive to the ring, has 



been interpreted as an X-ray irradiation effect on the HCN abundance (IKohno et al.l 12003 



20071 ). However, since the low- J lines only trace the cold and diffus e gas, higher- J lines may 



be better probes for the high ly excited nuclear environments (e.g., iMontero-Castano et al. 
20091 ; iMatsushita et al.ll2004l ). Here, we compare our first interferometric HCN(J = 3-2) and 
HCO + (J = 3-2) maps of NGC 1097, at an angular resolution of ~4", with the emission from 
the lower-J transitions. 



2. OBSERVATIONS AND DATA REDUCTION 



We observed NGC 1097 with the Submillimeter Arraj0 ( Ho et al. 2004 ). Two nights of 



observations were obtained in 2008 and 2010 with the compact configuration and the 345 
GHz receivers. The uv -coverages were ~7.5 - 61.5 k\ for both observations. The system 
temperatures (T sys ) were ~150-300 K and ~100-180 K, respectively, for the two sessions. 

The SMA correlator processes two IF sidebands separated by 10 GHz, with ~2 GHz 
bandwidth each in 2008. The bandwidth was upgraded to 4 GHz bandwidth for each IF after 
2009. The upper/lower sidebands were divided into slightly overlapping 24 and 48 chunks of 
104 MHz width in 2008 and 2010, respectively. In 2008, we observed the HCN(J = 3-2) line 
(rest frequency:265.886 GHz) in the LSB with 6 antennas. In 2010, we observed both the 
HCN(J = 3-2) and the HCO+(J = 3-2) (rest frequency:267.558 GHz) lines simultaneously 



1 The SMA is a joint pro ject between the Smithsonian Astrophysical Observatory and the Academia 
Sinica Institute of Astronomy and Astrophysics and is funded by the Smithsonian Institution and the 
Academia Sinica. 
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in the LSB with 8 antennas because of the availability of the 4 GHz bandwidth. 

We calibrated the SMA data with the MIR-IDL software package. The bandpass cali- 
brators are 3C84, Uranus and Ceres for the first track, and 3C454.3 and 3C84 for the second 
track. Callisto and Uranus were used for the absolute flux calibration for the first and sec- 
ond dataset, respectively. The accuracy of the absolute flux calibration is ~ 15% for both 
of the data. The gain calibrators are the quasars 0423-013 for the first track and 0137-245 
(phase) / 0403-360 (amplitude) for the second track, respectively. Mapping and analysis were 
done with the MIRIAD and the NRAO AIPS packages. We made the maps with natural 
weighting, and with the data binned to 20 km s _1 velocity resolution. The CLEAN processes 
were done in MIRIAD to remove the sidelobes. The synthesized beams are 3'.'7x2'. / 6 (PA = 
-10°) for the HCN(J = 3-2) maps and 4"4x2"6 (PA = -15°) for the HCO+(J = 3-2) maps. 
For the latter comparisons of both lines, we convolved the maps to 474x2 / . / 7 (PA = -14.6°). 
The la noise levels are llmJy beam -1 for both lines. 



Results 



3.1. High Excitation Dense Gas in the Circumnuclear Region of NGC 1097 



In FigureCDwe show the HCN(J = 3-2), HCO+(J = 3-2) and 12 CO(J = 3-2) flHsieh et all 
201 ll ) integrated intensity maps. The HCN(J = 3-2) and HCO + (J = 3-2) maps show a 
nuclear co ncentration and a molecular ring similar to the J = 1-0, J = 2-1, and the J = 3-2 



CO maps (IKohno et all 1200.1 120071 : iHsieh et all I2008L l201lf >. The nucleus is the brightest 
feature in all the maps, and the relatively weaker molecular ring coincides with the starburst 
ring. The ring consists of GMAs at the scale of giant molecular cloud associations (GMAs) 
of 200-300 pc. However, since both of the HCN and HCO + lines have higher critical density 
than CO, the HCN and HCO + GMAs are expected to be the denser regions of the GMAs 
mapped in the previous CO lines. We identify a total of 8 GMAs including the one at the 
nucleus (Tabled]). The GMAs are defined by the intensity peaks with visual inspection. We 
show the spectra of the nucleus in Figure [2j The spectra show asymmetric profiles, which 
are likely due to multiple underlying components. The blueshifted part is fainter than the 
redshifted part in all the lines. 

The nucleus is brighter in HCN than HCO + in both the low- and high-J lines. HCN(J 
= 3-2) has better contrast between the nucleus and the starburst ring than CO (J = 3- 
2) and HCO + (J = 3-2). We demonstrate this by explicitly comparing the contribution of 
the emission from the nucleus and the molecular ring in these lines. The HCN(J = 3-2) 
integrated flux for the nucleus (radius: 0"-6") and the ring (radius: 6"-13") are 34.4±3.4 Jy 
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km s 1 and 89.9±6.7 Jy km s l , respectively. Hence the nucleus contributes 30±3% of the 
HCN(J = 3-2) flux within 13". The HCO + (J = 3-2) integrated flux of the nucleus and the 
ring are 17.4±3.4 Jy km s _1 and 85.8±6.7 Jy km s _1 , respectively. Therefore, the nucleus 
contributes 17±3% of the HCO+(J = 3-2) flux within 13". The integrated 12 CO(J = 3-2) 
flux of the nucleus and ring are 1021 and 3778 Jy km s -1 . Hence the nucleus contributes 
21±0.2% of the total 12 CO(J = 3-2) flux within 13". The HCN(J = 3-2) is therefore 
significantly brighter in the nuclear region as compared with the molecular ring. 



3.2. Intensity Ratios of the Dense Gas in the Circumnuclear Region 



We compare the HCN(J =3-2), HCO+(J = 3-2), and 12 CO(J = 3-2) flHsieh et al.ll201lh 
at the angular resolution of 4"4x2"7. For the comparison of the J = 3-2 and J = 1-0 lines, 
we limited the shortest baselines to 7.5 kX for all transitions, in order to match the resolution 
of the J = 1-0 lines at 7"9x3"2. HCN(J =1-0), HCO+(J = 1-0), and 12 CO(J = 1-0) lines 



Kohno et al. 


2007, 


2003. 


2011) 



In TableH the averaged intensity ratios of HCN(J = 3-2)/HCO + (J = 3-2) are 1.11±0.36 
and 2.13±0.34 in the ring and the nucleus, respective ly. These values are similar to the 
single dish observations compiled by iKrips et al.l ( 120081 ). where the ratio is nearly unity in 
the starburst galaxies and larger than unity in the galaxies which host AGNs. The average 
intensity ratios of HCN(J = 3-2)/ 12 CO(J = 3-2) are 0.05±0.01 and 0.07±0.01 in the ring 
and the nucleus, respectively. In Figure [3] we show the correlation of the fluxes of HCN(J = 
3-2), HCO+(J = 3-2), and 12 CO(J = 3-2) lines. While there is a general correlation, the 
dispersion of these ratios within the molecular ring is about a factor of 2. We calculated the 
linear Pearson correlation coefficients of the plots in Figure |3J The correlation coefficients 
of all data are 0.86, 0.69, and 0.63 for Figure Eta), (b), and (c) respectively. The correlation 
coefficients of the GMAs in the starburst ring are 0.74, 0.44, and 0.21 for Figure^a), (b), and 
(c). The higher coefficients represent stronger linear dependence of the two quantities. This 
variation in the molecular ring is responsible for the apparent differences between the GMAs 
in the ring. Such differences are marginally significant but difficult to interpret. Chemistry, 
excitation, and the contributions from diffuse gas to the CO(J = 3-2), may all be important. 

We compare the ratios of J = 3-2 and J = 1-0 lines in low resolution (7"9x3"2) in 
Table[2 iKohno et all toOlh reported that the lower-J line ratios of HCN(J = 1-0)/ 12 CO(J = 
1-0) is -0.34, and the HCN(J = 1-0) /HCO+(J = 1-0) is ~2 in the nucleus. The HCN/HCO+ 
seems to be similar in lower- and higher-J lines, while the HCN/ 12 CO decreases toward the 
higher-J lines in the nucleus (0.04±0.01). The ring shows similar trends in Table [2j In 
both nucleus and ring, the HCO + /CO ratios decrease in J = 3-2. On the other hand, the 
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HCN(J = 3-2)/HCN(J = 1-0) and HCO+(J = 3-2)/HCO+(J = 1-0) ratios does not show 
significa nt difference i n the ring and the nucleus. This result is consistent with the survey 
done by iKrips et al.l ( 120081 ). The decreasing HCN/CO ratio s toward higher- J lines were 
repor ted in other galaxies hos t AGN, for i nstance, NGC 1068 (IKrips et al.l 120081 1201 if ), and 
M51 JMatsushita et al.l Il998l . l2004i wj). We will discuss the physical conditions of the 
molecular gas in the Sect. HI 

We derive the line ratios in LTE conditions to roughly look into the effects of temper- 
atures, opacities and abundance ratios on the HCN and HCO + lines. A more sophisticated 
radiative transfer model will be important once we have more lines and higher angular res- 
olution. In our case, we are looking at the large scale properties of the molecular gas, and 
a simple LTE analysis is adequate since all the information are averaged by the coarse res- 
olution. We define the brightness temperature T B = /[J„(T ex ) — J u (T hg )](l — e _T ), where 
Tb is the brightness temperature, / is the filling factor, J v = (hv/k)(l/e hu / kT ~ 1 ), T ex is the 
excitation temperature, Tb g is the background temperature, and r is the opacity. In our 
estimates we assume the filling factors are the same for both lines. For the same molecular 
species at J = 3-2 and J = 1-0, the estimates are: 



1. Both lines are optically thick (r ^> 1), in which case the integrated intensity line ratio 
can be written as 

jr>thick _ Wv32 C^ex) ~ ^32 (^bgjj /-|\ 
[Jv 10 (T cx ) — Ju 10 {Tbg)] 

where the background temperature T bg is 2.7K. In LTE, the -R32/10 approaches to unity 
with rising temperatures (see Figure and Figure [5^). 

2. Both lines are optically thin (r <C 1) and the integrated intensity line ratio can be 
written as 

ithin [Jv32 C^ex) ~ ^32(^bg)]T32 



R 



32/10 



[J Vl0 (T ex ) - Jv 10 (T bg )]rio' 



and 



732 



,^32 



hv io 1 



3 -/W32/fcT ex 



— = 3 — e kT <* 
T10 ^10 1 



(2) 



(3) 



e -/i^io/fcTex ' 

which in LTE depends on the optical depth r at each frequency. r 32 and r 10 are the 
opacity of J = 3-2 and J = 1-0 lines. In Figures and Eh, the ratio increases rapidly 
with T ex . 



3. We considered J = 3-2 line is optically thin and J = 1-0 line is optically thick 



R. 



thin/thick 
32/10 



[J i/32 (T ex ) ^; i2 (^~bg)] 
WviO (-^ex) — ^fio(-^bg)] 



T32- 



(4) 
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In Figure [lb and Figure [5b, we illustrate the dependence of the ratios on temperature 
and opacity, where the ratios of the same molecules are independent of the abundance. 
In both figures, the ratios seem to be more sensitive to the temperatures when T ex < 
20K, but are dependent on the opacity of HCN(J = 3-2) and HCO+(J = 3-2) at T cx > 
20K. This dependence of temperature/opacity also varies with ratios. 

The ratios of 12 CO(J = 3-2)/ 12 CO(J = 1-0) lines are shown in Tabled We also derive the 
LTE results of the 12 CO(J = 3-2)/ 12 CO(J = 1-0) ratios in Figure El The first two cases 
in Figure Eh are for when both lines are optically thick and optically thin. The third case 
(Figure [6b) is assuming optically thin 12 CO(J = 1-0) and optically thick 12 CO(J = 3-2) 
lines. 

For the ratios of the different species of HCN(J = 3-2) and HCO + (J = 3-2), we illustrate 
the two cases when both lines are optically thick and thin: 

1. Both lines are optically thick (r» 1), in which case the integrated intensity line ratio 
can be written as 

Dthick _ [JHCN{T ex ) - Jgcv(T b g)] 

- fl HCN/HCO+ — 77 7^f~\ 7 ~Fr w ' v > 

l'JHCO+{ I c X ) — •JHCO+\ 1 hg)\ 

The ratios quickly saturate to nearly unity (see Figure Uh- and Figure Uk>) in the case 
when both lines are thick. 

2. Both lines are optically thin (r <C 1) and the integrated intensity line ratio can be 
written as 

^thin + = [JhCn(T cx ) - JHCN(T hg )}T HCN 

HCN/HCO [JHCO+(Tex) — JHCO+{Tbg)} T HCO+'' 



and 



9 ~ ~ 6hB HCN32 . n rp 

THCN _ ^HCO + 32 ^HCN -B HC N32 e ^ A HC N32 1 - e - h "HCN/kT e 



(7) 

^HCO+ ^HCN32 ^HCO+ ^HCO + 32 e ^t^ 32 ^HCO+32 1 - e-' w HCO+ / tT - ' 

which in LTE depends on the optical depth r at each frequency. t H cn and r H co+ 
represent the opacity of HCN(J = 3-2) and HCO + (J = 3-2) lines. The Shcn32 and 
-E?hco+32 are the rotational constants of HCN and HCO + . Ahcn32 and A H co+32 are the 
Einstein A coefficients of HCN(J = 3-2) and HCO+(J = 3-2). A^hcn and A^ HC o+ are 
the total column density of HCN and HCO + molecules. Therefore the free parameters 
are the column density ratio, which is proportional to the abundance ratio of the HCN 
and HCO + molecules if their size scales and the excitation temperature T ex are similar. 
In Figure [7b, we show the HCN(J = 3-2)/HCO + (J = 3-2) ratios as a function of 
T ex and HCN/HCO + abundance ratio. The line ratios are nearly insensitive to the 
temperatures, but are dependent on the abundance ratios. 
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3.3. Correlation of the molecular line intensity and 24/mi continuum 

In Figure [8] we show the comparison of the molecular gas and the archival Spitzer 24//m 
continuum emission. The PSF of the 24/zm map is ~6", and we smoothed our molecular 
line maps to match this resolution. In Figure [9] we show the relation of the molecular line 
intensity with the 24/xm continuum emission for the GMAs. In Figure E^a), the HCN(J = 
3-2) flux of the ring seems to be tightly correlated to the 24/zm continuum emission, while 
the nuclear emission is offset from this correlation. The relation is less tight for the lower 
critical density lines of HCO + (J = 3-2) and 12 CO(J = 3-2). However, for these lines, the 
nuclear emission is closer to the correlation of the ring. We calculated the linear Pearson 
correlation coefficients of the plots in Figure M The correlation coefficients of all data are 
-0.23, 0.26, and -0.10 for Figure 12(a), (b), and (c) respectively. The correlation coefficients 
of the GMAs in the starburst ring are 0.86, 0.51, and 0.36 for Figure [9(a), (b), and (c). In 
this case, including the nucleus will decrease the linear dependence of the two quantities. 



DISCUSSIONS 



The nucleus of NGC 1097 seems to have environments of higher excitation than the ring. 
Our spatially resolved HCN(J = 3-2) and HCO + (J = 3-2) maps show that the nucleus is 
brighter than the molecular ring. Similar results have also been shown in o ther interferomet- 



ric studies o f mul t i-J HCN/HCO + e mission such as M51 and NGC 1068 ( [Matsushita et al. 



19981 . 12004 . 120071 ; iKrips et al.l 1201 ll ). The ratio of upper transition as compared to lower 
transitions, on the other hand, does not show a sharp enhancement in the nuclear region. 
The nature of the higher excitation in the nucleus can thus be deduced as following discus- 
sions. The necessity o f analysis wit h mult iple transitions has already b een indicated from 
the previous studies by iKohno et al.l ( 120031 ) and the single dish survey by IKrips et al.l (120081 ) 
toward nearby active/starburst galaxies. In this paper, we discuss our results with standard 
LTE analysis. 



4.1. Physical Properties of the J = 3—2 lines 

In Sect. 13.21 we derived the LTE estimates of the physical parameters for the HCN and 
HCO+ molecules. In Figures Si andEfc, if the HCN(J = 3-2), HCN(J = 1-0), HCO+(J = 3- 
2), and HCO + (J = 1-0) lines are optically thick, which suggests cold/thick dense molecular 
gas in both the nucleus and the ring with T ex ~ 5K. This is possible since our spatial 
resolution is still quite coarse at 300 pc, so that we are likely sampling extended and colder 
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material. However, sin ce the optically thick molecular gas has typical temperature of >10K 



flDishoeck et all 11993 ). while the observed HCN(J = 3-2)/HCN(J = 1-0) and HCO+(J = 
3-2)/HCO + (J = 1-0) ratios are lower than LTE estimates of ~ 0.5 at 10K. The inconsistence 
may be either the four lines are not all optically thick, or the different filling factors of those 
transitions. 

In Figure [TJi, the HCN 32 /HCO;j 2 ratio is nearly unity if both lines are thick. This is 
consistent with the average ratio of the ring, but inconsistent for the nucleus. Therefore, both 
HCN(J = 3-2) and HCO + (J = 3-2) being optically thick is less likely for the nucleus. This 
case can not be ruled out since we do not have a promising optically thin line to determine 
the physical parameters. Alternatively, that the J = 3-2 lines are weaker than the J = 1-0 
lines, may be because the J = 3-2 lines are thin and J = 1-0 lines are thick. As it is shown in 
Figures Sb and [5b, the line ratios depend on the J = 3-2 opacity and T ex . Note that the J = 
3-2 to J = 1-0 ratios for the nucleus and the ring, are similar within the uncertainties. The 
LTE estimates suggest that the J = 3-2 lines may be colder with higher opacity or warmer 
with lower opacity. 

In the case where both HCN(J = 3-2) and HCO + (J = 3-2) are optically thin, we could 
estimate the abundance ratio of HCN/HCO+ in Figure [7b. The [HCN]/[HCO+] ratio is 
~ 2.5 — 4.5 and ~ 1 — 3 for the nucleus and the ring. The optically thin scenarios seem to be 
supported by Figure [9k,b, which show tighter correlations with the 24/zm emission at GMAs 
scale. The tighter relations of HCN(J = 3-2), and HCO + (J = 3-2) with 24/im emissions 
suggest these two lines trace the star forming regions better than CO (J = 3-2) line. If the 
HCN(J = 3-2), HCO + (J = 3-2), and 24/zm emissions are optically thin, their fluxes will be 
proportional to the column density of molecular gas and dust. Their good correlations may 
reflect that these are reliable tracers of the mass. On the other hand, the abundant molecule 
of CO is likely to remain optically thick in the J = 3-2 line, which leads to a poor relation 
in Figure [St. We note that the enhancement of the HCN(J = 3-2) emission in the nucleus is 
not simple. The HCN(J = 3-2) has a factor of two excess flux as compared to the HCO + (J 
= 3-2). This excess is likely due to the different abundance ratio of the HCN and HCO + in 
the ring and the nucleus as it is shown in Figure [7b. The abundance ratio would depend on 
the different physical environments, and we will discuss the effects of the chemistry in the 
following section. 

On the contrary to the low HCN(J = 3-2)/HCO+(J = 3-2) ratios, the high 12 CO(J = 
3-2)/ 12 CO(J = 1-0) ratios in the circumnuclear region of NGC 1097 seem to suggest that 
the 12 CO(J = 3-2) line is optically thick, and the 12 CO(J = 1- 0) is optically thin (Figure ED. 



The similar hig h ratios >2 were also r eported in NGC 1068 (IKrips et al.ll2011t iTsai et al. 



201lh and M51 (Matsushita et al.ll2004f ). In Figure El for the case of 12 CO(J = 3-2)/ 12 CO(J 
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= 1-0) ratios larger than unity, the higher T ex or lower opacity of 12 CO(J = 1-0) line lead to 
the higher 12 CO(J = 3-2)/ 12 CO(J = 1-0) ratios. Given the similar 12 CO(J = 1-0) opacity 
of the nucleus and the molecular ring, the higher 12 CO(J = 3-2)/ 12 CO(J = 1-0) ratio of the 
nucleus suggests it has higher temperature than the ring, or the nucleus has lower 12 CO (J 
= 1-0) opacity than the ring. These factors might also be consistent with the decreasing 



HCN /CO ratios toward high- J lines reported in the previous studies by iKrips et al.l ( 12008 



201ll ). Note that the LTE analysis gives the simplified results since the ratios also depend 
on other physical parameters such as gas density. Unfortunately, the LTE analysis can not 
provide straightforward density solution. This is because we need additional information 
such as fractional abundance and source size. However, in principle we expect in the case of 
collisional excitation, higher density conditions might increases the higher-J population of 
molecules. 

In terms of excitation, the observed low values for the ratios of HCN(J = 3-2)/HCN(J 
= 1-0) and HCO + (J = 3-2)/HCO+(J = 1-0) also suggest that the J = 3-2 transitions 
might be subthermal. The case of subther mal excitation for the higher- J HCN/HCO" 1 " line s 



has also been found in other galaxies (e.g., iNguyen-Q-Rieu et al.lll992l ; iRiechers et al.ll201ll ). 



The critical density is defined when the molecules are in thermal equilibrium via collisions. 
Densities have to be higher in order to drive the collisional time scale to be significantly 
shorter than the spontaneous decay time scale as defined by the Einstein A coefficient, in 
order to populate the upper energy levels. Since the Einstein A coefficient scales as z/ 3 , the 
rapid spontaneous decay rate will eventually overwhelm the collision rate for a sufficiently 
high energy state. In such case, the subthermally excited lines will set an upper limit on 
the densities as compared to the critical density of the particular transitions. Based on 
the bright HCN/HCO + J = 3-2 emission, and their low intensity ratio as compared to 
the lower transitions, the low density scenario with subthermal excitation, is also possible. 
However, the high 12 CO(J = 3-2)/ 12 CO(J = 1-0) ratios of the circumnuclear region in 
NGC 1097 might suggest high gas density, which is in conflict to the above scenario. A 
careful excitation analysis cooperating with multi-J molecular lines are essential to derive 
the gas density, excitation temperature, and kinetic temperature to examine the possibility of 



subth e rmal conditions of J = 3-2 HCN, HCO + lines, such as the works done by lKrips et al. 



(|2008l . 1201 if ) 



4.2. Column Density of the HCN and HCO + molecules 



We can calculate the column density of HCN and HCO + by assuming LTE conditions 
and correcting for the optical depth, 
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3k e^roM+mn* + hBrot/3k) Ty r 
mo1 8^B rot ^ J t + 1 (hu/k)e~^/kT ex !_ e -r„ J W 

where A mo i is the column density of specific molecule, k and h are the Boltzmann and Planck 
constants, B TOt is the rotational constant of the molecule, fi is the electric dipole moment of 
the molecule, J/ is the quantum number of the lower energy level from J + 1 — > J, T ex is the 
excitation temperature which is constant under LTE assumption for all transitions, r„ is the 
opacity, and dv is the velocity width. The column density of H 2 can then be derived from 
the abundance ratio of [Molecule] /[H 2 ]. Alternatively, if we derive the column density of H 2 
from other observations, we can estimate the abundances of HCN and HCO + . We calculate 
B TOt of HCN and HCO+ to be 44.31 and 44.59 GHz. The values of pi for HCN and HCO+ 
are 2.98 xl0~ 18 and 3.92 xl0~ 18 e.s.u. cm. We then calculate A H cn and iV H co+ by 

iV H cN = 1-108 x 1^1^^-^— J T b (HCN(J = 3 - 2))dv, (9) 
AW = 6.325 x 10 9 ( ^ 25 + 68 ° 5 ^\_ r ^ / T b {HCO + {J = 3 - 2))dv, (10) 



We estimate the cases of the nucleus and the ring separately with different temperatures 
and opacities. As for the nucleus, we adopt the warm and thin c o nditio ns based on the 



previous discussions. In terms of temperatures, ISandstrom et al.l (120101 ) fitted the SED 
of the nucleus of NGC 1097 for the spatial scale of 600 pc, and found that the thermal 
blackbody spectrum peaks at ~ lOO/xm, corresponding to a dust temperature of Td ~ 30K. 
However, an additional peak seems to be located at ~10/im, where corresponds to T d of 
~ 300K. We thus constrain the range of the A H cn within this range. Since the opacities 
are nearly constant from 30 to 300 K (Figure Hb), the corresponding opacities of HCN(J 
= 3-2) and HCO+(J = 3-2) are ~0.2 in this range. Then A HC n is ~2xl0 13 to lxlO 14 
cm -2 , and Ahco+ is ~6xl0 12 to 3xl0 13 cm" 2 . The column density of H 2 in the nucleus, as 
derived from CO (J = 3-2) (Tabled]) with a conversion factor of 1.8xl0 20 (K km s -1 ) -1 cm -2 



(IDame et al.ll200ll ). is 3xl0 22 cm 2 . If we adopt this value for the total H 2 column density, 
we can estimate the fractional abundance of HCN and HCO + . We find Ahcn/Ah 2 is 7x 10~ 10 
to 3xl0" 9 , and A HC o+/A H2 is 2x10 10 to 9x10 10 . The conve rsion factor in the galactic 



center was suggested to be lower than the value of galactic disk. iMauersberger et al.l (119961 ) 
suggests that the conversion factor is 0.3 x 10 20 (K km s -1 ) -1 cm -2 in the nuclear starburst 
environments. Hence the abundances mentioned above would be larger by a factor of 6. 
The key uncertainty here is whether the CO, HCN, and HCO + , are spatially coextensive. 
The derived column density ratios are the mean values averaged over the volume of the CO 
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emission. Another route to estimate the [HCN]/[H 2 ] is to adopt the a s deri ved from the 
X-ray absorption. The Nn 2 of the nucleus as derived by iNemmen et al.l ( 120061 ) is 4.6 xlO 20 



cm 2 , which increases the [HCN]/[H 2 ] and [HCO + ]/[H 2 ] by a factor of ~60. 

As for the molecular star-forming ring, using the same analysis, with a temperature 
of 30 K, and opacities of 0.2, we derive Ah 2 for the ring-GMAs of 3xl0 22 cm -2 , Ahcn of 
~7x 10 12 cm" 2 , and iV HC o+ of ~ 4 x 10 12 cm" 2 . The [HCN]/[H 2 ] and [HCO+]/[H 2 ] are 
2xl0~ 10 and lxlO -10 , respectively. 

CO is usually optically thick because of its large abundance and low critical density. The 
concern is that we may sample the cold gas on the surface of the molecular clouds with CO, 
while the HCN and HCO + sample the inner region of the clouds. The H 2 column density 
as derived from the X-ray absorption may be more appropriate if the X-ray traces the more 
compact and denser regions. We thus favor the larger fractional abundance as estimated 
above via the X-ray absorption. 



4.3. Chemical Abundances in the Nucleus 



We can check if these chemical abu ndances are consis t ent w ith the ionization chemistry 
as suggested in the classical studies of iLepp &: Dalgarnol (119961 ) . These authors predicted 
the fractional abundance of molecules for various ionizing densities (see their Fig. 2, 3). The 
ionizing flux densities was defined by £/nn 2 (erg s -1 cm -3 ) in their paper, where £ and nn 2 are 
the ionization rate and the volume number density of the molecular hydrogen. If the ionizing 
flux is dominated by the X-ray emi ssion in the nucleus, i ts luminosity with i n the energy range 
of 2-10 KeV is 4.4 xlO 40 erg s" 1 jNemmen et alJbood ). llVlalonev et aD Jl996h derived the 



X-ray ionization rate in the low-ionization limit as £ ~ 1.4 x 10- n Lur2 2 Nz2 (s _1 ), where 
Lu, r 2i A 22 are the X-ray luminosity in unit of 10 44 erg s _1 , the distance to the X-ray source 
in the unit of 100 pc, and the H 2 column density attenuating the X-ray flux in the unit of 
10 22 cm -2 . We assume the ionization radiation has an isotropic distribution, and we adopt 
the -/Vh 2 as derived from CO of 3xl0 22 cm -2 at the radius of l'/6 = 126 pc (the average half 
width of the synthesized beam). This gives £ of 1.3 xl0~ 15 s _1 . We have also estimated that 



fiH 2 of the nucleus is 10 4±1 cm 3 (IHsieh et al.l 120081 ). Hence log£/nn 2 is -17.9 to -19.9. For 



the other case of the Ah 2 estimated from the X-ray absorption, the corresponding £ is 10~ 13 
s _1 and the log£/n H2 is -15.9 to -17.9. The HCN fractional abundance predicted by the 
model corresponds to the range estimated by CO and X-ray absorption are -8.7 to -7.4 and 
<-7.8, respectively. Thus, the HCN abundance predicted by the X-ray ionization chemistry 
is consistent with both cases we estimated above. However, if it favors the latter case, a 
larger H 2 number density is probably required to generate the higher HCN abundance. We 
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also find that the HCO + abundance is consistent with the X-ray ionization chemistry for 
the case of Nn 2 as derived from X-ray absorption. We conclude from our J = 3-2 high 
resolution studies, that the X-ray ionization chemistry can explain the physical conditions 
in the nucleus. 



5. SUMMARY 

1. We detected the J = 3-2 HCN, HCO+ lines in the circumnuclear region of NGC 1097. 
The HCN and HCO + lines reveal a bright nucleus as compared to the molecular ring. 
The HCN(J = 3-2) lines contribute ~30% of the total flux to the nucleus, while it 
is only -20% for the HCO+(J = 3-2) and CO(J =3-2) lines. The HCN(J = 3- 
2)/HCO + (J = 3-2) has a value of 2 in the nucleus and unity in the ring. The J = 3-2 
to J = 1-0 ratios of both lines are as low as ~0.2 for both the nucleus and the ring. 
Our results generally agree with the previous single dish surveys. 

2. Our standard LTE analysis suggests that the HCN (J = 3-2) and HCO + (J = 3-2) lines 
might be optically thin in the nucleus. Both lines can be optically thin or thick in 
the molecular ring. The high-J lines with larger critical density, however, might also 
be subthermally excited. The tight correlation of the flux of HCN and HCO + and 
24/zm emission suggest these two lines trace the same origins of star formation in the 
ring better than CO, while other mechanisms such as AGN activities contribute to the 
HCN enhancements in the nucleus. 

3. The simple estimates of the self-consistency of the X-ray ionization chemistry is made 
with our high-J observations. The results suggest that the ionization chemistry model 
can explain the excess of the HCN abundance in the nucleus. 
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Table 1. Ratios at high angular resolution 



GMAs 


Positions 


HCN(J = 3-2) 
(1) 

(K km s- 1 ) 


HCO+(J = 3-2) 
(2) 

(K km s- 1 ) 


12 CO(J = 3-2) 
(3) 

(K km s- 1 ) 


HCN32/HCO+32 
(4) 


HCN 3 2/ 12 C0 3 2 
(5) 


HCO+32/ 12 C0 32 
(6) 


1 


3'.'0,9?5 


10.8 


8.1 


189.0 


1.33±0.34 


0.05±0.01 


0.05±0.01 


2 


-5'.' 0,7'.' 5 


7.6 


5.2 


120.9 


1.47±0.57 


0.06±0.01 


0.05±0.01 


3 


-8'.'5,-3'.'0 


9.5 


8.5 


273.4 


1.11±0.29 


0.03±0.01 


0.04±0.01 


4 


-5('0,-8'.'0 


7.9 


7.5 


159.5 


1.06±0.32 


0.04±0.01 


0.04±0.01 


5 


l'.'0,-9'.'0 


6.4 


4.8 


110.4 


1.36±0.59 


0.06±0.02 


0.06±0.02 


6 


8'.' 5, -6'.' 5 


6.0 


9.7 


121.9 


0.62±0.20 


0.04±0.01 


0.06±0.01 


7 


8"5,4'.'0 


8.1 


9.7 


175.3 


0.84±0.22 


0.04±0.01 


0.05±0.01 


8 


0'.'5,0'.'0 


24.7 


11.6 


336.1 


2.13±0.34 


0.07±0.01 


0.04±0.01 



Note. — The flux density of the peaks are measured within one synthesized beam (4"4x2'.'7). The (<5R.A.,<5Decl.) position 
of the peaks relative to the phase center. The uncertainties of the HCN(J = 3-2), HCO+(J = 3-2), and CO(J = 3-2) are 1.7 
K km s , 1.7 K km s~ 1 , and 3.6 K km s , respectively. 



Table 2. Ratios at low angular resolution 



Ratios Nucleus Ring 

(1) (2) 



HCNio/HCO+io 2.22±0.50 1.67±0.83 

HCNio/COio 0.59±0.13 0.35±0.14 

HCO+io/COio 0.27±0.07 0.21±0.11 

HCN32/HCO+32 2.01±0.41 1.07±0.48 

HCN32/CO32 0.04±0.01 0.02±0.01 

HCO+32/CO32 0.02±0.01 0.02±0.01 

HCN32/HCN10 0.15±0.02 0.12±0.05 

HCO+32/HCO+10 0.17±0.05 0.19±0.10 

12 CO 3 2/ 12 COi 2.33±0.44 1.84±0.56 



Note. — In this table, we show the 
low resolution (7"9x3'.'2) intensity ratios 
of the nucleus and ring. The ratio of the 
nucleus is measured at the intensity peak. 
We measured the azimuthally averaged 
intensity ratios of the ring within radius 
of 6"-13". 
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Table 3. Comparisons of flux of 24/zm and molecular lines 



Peaks 


HCN(J = 3-2) 
(K km s- 1 ) 


HCO+(J = 3-2) 
(K km s- 1 ) 


12 CO(J = 3-2) 
(K km s- 1 ) 


24/im 
(mjy arcsec -2 ) 


1 


5.0±0.8 


5.5±0.8 


100.2±1.8 


114.49±1.64 


2 


4.1±0.8 


3.6±0.8 


72.2±1.8 


113.82±1.64 


3 


4.3±0.8 


4.4±0.8 


146.7±1.8 


102.97±1.56 


4 


3.9±0.8 


4.7±0.8 


100.4±1.8 


82.95±1.40 


5 


4.0±0.8 


3.1±0.8 


75.7±1.8 


77.99±1.35 


6 


3.2±0.8 


2.9±0.8 


72.4±1.8 


56.66±1.15 


7 


3.9±0.8 


5.4±0.8 


99.0±1.8 


85.99±1.42 


8 


12.3±0.8 


5.6±0.8 


190.1±1.8 


68.91±1.27 



Note. — We show the low resolution (6") intensity of the nucleus and 
ring, as well as the Spitzer 24/im continuum emission. The unit of the flux 
density of 24/im continuum emission is mjy arcsec - 2 (10 -3 erg s~ 1 cm -1 
arcsec -2 Hz -1 ) integrated over 6". 
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10 -10 10 -10 10 -10 



RA offset (arcsec; J2000) RA offset (arcsec; J2000) RA offset (arcsec; J2000) 

Fig. 1. — Left: HCN(J = 3-2) integrated intensity map, the contours levels are, 2, 4, 6, 8, 
10, 12, 14cr, where la is 1.2 Jy beam -1 km s -1 . Middle: HC0 + (J = 3-2) integrated intensity 
map, the contours levels are, 2, 4, Qa, where la is 1.2 Jy beam -1 km s -1 . Right: 12 CO(J = 
3-2) integrated intensity map, the contours levels are, 2, 4, 8, 16, 32, 64a, where la is 4.5 
Jy beam -1 km s -1 . The synthesized beam of the three images are 4 / .'4x2'. / 7 (PA = -14.6°). 




Fig. 2.— (a) HCN(J = 3-2) spectrum of the nucleus (solid-line). We subtract the HCN(J = 
3-2) spectrum by CO spectrum (Fig. 2c) as template to check whether the residuals (dotted- 
line) are significant. The CO spectrum is multiplied by 0.05 for the subtraction, (b) HCO + (J 
= 3-2) spectrum of the nucleus (solid-line). We subtract the HCO + (J = 3-2) spectrum by 
CO spectrum as template to check whether the residuals (dotted-line) are significant. The 
CO spectrum is multiplied by 0.03 for the subtraction, (c) CO(J = 3-2) spectrum of the 
nucleus. The synthesized beam of all the spectra are 4'.'4x2 / . / 7 (PA = -14.6°) 
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Fig. 3.— (a) The correlation of HCN(J = 3-2) and 12 CO(J = 3-2) line intensity. The 
fluxes of the GMAs in the ring are shown as squares, and star represents the data of the 
nucleus. The units are K km s -1 . The ±ler uncertainties are overlaid on the symbols with 
vertical/horizontal bars. The linear Pearson coefficients are shown in the plots. r(AGN+SB) 
are the coefficients of all data, and r(SB Ring) are derived for the starburst ring GMAs only, 
(b) The correlation of HCO + (J = 3-2) and 12 CO(J = 3-2) line intensity, (c) The correlation 
of HCN(J = 3-2) and HCO + (J = 3-2). The symbols are the same as (a). The data we show 
here are from resolution of 4 / /4x2 / .'7. 
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Fig. 4. — (a) The HCN(J = 3-2)/HCN(J = 1-0) intensity ratio as a function of excitation 
temperature in LTE conditions. The solid line is for when both lines are optically thick. The 
dotted line is for when both lines are optically thin, (b) The HCN(J = 3-2)/HCN(J = 1-0) 
intensity ratio as a function of excitation temperature and thcn32 m LTE conditions. In this 
case, the HCN(J = 3-2) is thin and HCN(J = 1-0) is thick. The ratios are shown as contours 
from 0.05 to 0.5 in steps of 0.05. (c) The HCN( J = 3-2) /HCN( J = 1-0) ratios of the nucleus 
0.15±0.02 shown as a function of the T ex and r H cN 32 . (d) The HCN(J = 3-2)/HCN(J = 1-0) 
ratios of the ring 0.12±0.05 shown as a function of the T ex and t H cn 32 - 
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Fig. 5. — (a) The HCO + (J = 3-2)/HCO + (J = 1-0) intensity ratio as a function of excitation 
temperature in LTE conditions. The solid line is for when both lines are optically thick. The 
dotted line is for when both lines are optically thin, (b) The HCO + (J = 3-2)/HCO + (J = 
1-0) intensity ratio as a function of excitation temperature and t hco + in LTE conditions. In 
this case, the HCO + (J = 3-2) is thin and HCO + (J = 1-0) is thick. The ratios are shown as 
contours from 0.05 to 0.5 in steps of 0.05. (c) The HCO + (J = 3-2)/HCO + (J = 1-0) ratios 
of the nucleus 0.17±0.05 shown as a function of the T ex and t"hco + ■ 

(d) The HCO+(J = 

3-2)/HCO + (J = 1-0) ratios of the ring 0.19±0.10 shown as a function of the T ex and t hco + . 
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Fig. 6. — (a) The 12 CO(J = 3-2)/ 12 CO(J = 1-0) intensity ratio as a function of excitation 
temperature in LTE conditions. The solid line is for when both lines are optically thick. 
The dotted line is for when both lines are optically thin, (b) The 12 CO(J = 3-2)/ 12 CO(J 
= 1-0) intensity ratio as a function of excitation temperature and rcoio i n LTE conditions. 
In this case, the 12 CO(J = 1-0) line is thin and 12 CO(J = 3-2) line is thick. The ratios are 
shown as contours higher than unity, from 1.5, 3, and 4.5 to illustrate the dependence of 
the opacity and T ex . (c) The 12 CO(J = 3-2)/ 12 CO(J = 1-0) ratios of the nucleus 2.33±0.44 
shown as a function of the T ex and t COio . (d) The 12 CO(J = 3-2)/ 12 CO(J = 1-0) ratios of 
the ring 1.84±0.56 shown as a function of the T ex and t COio . 
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Fig. 7. — (a) The HCN(J = 3-2)/HCO + (J = 3-2) intensity ratios as a function of excitation 
temperature in LTE conditions. In this case both lines are optically thick. The ratios 
are nearly unity, (b) The HCN(J = 3-2)/HCO + (J = 3-2) intensity ratios (contours) as a 
function of excitation temperature and abundance ratios of HCN/HCO + in LTE conditions. 
The ratios of the contours are labeled. In this case, both lines are optically thin. 
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Fig. 8. — Left: The HCN(J = 3-2) map is overlaid on the 24//m continuum map. Middle: 
The HCO + (J = 3-2) map is overlaid on the 24/im continuum map. Right: The CO(J = 
3-2) map is overlaid on the 24/im continuum map. The resolutions of the HCN(J = 3-2) 
and HCO + (J = 3-2) are matched to that of the PSF of the 24/zm continuum map (6"). The 
contour starts from 4.50 Jy Beam -1 km s , and in steps of 1.50 Jy Beam -1 km s -1 . 
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Fig. 9. — (a) The correlation of HCN(J = 3-2) integrated intensity and 24//m continuum 
flux. The symbols are the same as Figure [3j The drier uncertainties are overlaid on the 
symbols with vertical/horizontal bars. The linear Pearson coefficients are shown in the plots. 
r(AGN+SB) are the coefficients of all data, and r(SB Ring) are derived for the starburst 
ring GMAs only, (b) The correlation of HCO + (J = 3-2) integrated intensity and 24/zm 
continuum flux. The symbols are the same as Figure [3j (c) The correlation of 12 CO(J = 
3-2) integrated intensity and 24/zm continuum flux. The symbols are the same as Figure |3j 
The flux are measured from low resolution maps (6"). 



